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Abstract 

The forces acting on and the energies of solitons governed by the 
nonlinear Schrodinger equation in finite, planar waveguides with peri- 
odic and with homogeneous Dirichlet, Neumann and Robin boundary 
conditions are determined by means of a quantum analogy. It is shown 
that these densities have S-shape profiles and increase as the hardness 
of the boundary conditions is increased. 
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1 Introduction 



Soliton propagation has received a great deal of attention in recent years 
because of its possible applications in distorsionless signal transmission in 
ultra-high speed and long-distance transoceanic telecommunications by op- 
tical fibers [T|, soliton lasers, all-optical soliton switches [2], nonlinear pla- 
nar waveguides [3], etc. However, most of the theoretical studies on one- 
dimensional solitons have been concerned with the initial- value problem, i.e., 
infinite intervals in space, or with boundary-generated solitons, i.e., semi- 
infinite intervals, and have been based on the inverse scattering transform [I] , 
perturbation methods or the infinite sine Fourier transform [5]. 

Since nonlinear transmission lines, optical fibers, etc., have a finite length, 
the infinite and semi-infinite spatial domains used in most previous studies 
of nonlinear wave propagation are idealizations. In this letter, the nonlin- 
ear Schrodinger equation in finite intervals with periodic and homogeneous 
Dirichlet, Neumann and Robin boundary conditions is studied numerically 
in order to determine the energy of and the forces on solitons prior to, during 
and after their interaction with the boundaries. The energy of and the forces 
on solitons are evaluated by means of a quantum analogy between wave prop- 
agation phenomena governed by the nonlinear Schrodinger equation and the 
motion of a particle in a nonlinear potential well. Bian and Chan [3J have 
previously studied the nonlinear and diffraction forces on solitons governed 
by the initial-value problem of the nonlinear Schrodinger equation. These 
authors, however, used a filter and their force and energy densities are not 
consistent with a strict quantum mechanics analogy. 
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2 One-Dimensional, Nonlinear Schrodinger Equa- 
tion 



The envelope of the electric field in a two-dimensional medium with a quadratic 
nonlinear refractive index, n=n +n 2 \E'\ 2 , i.e., subject to the Kerr effect, with 
neither losses nor other higher order effects, can be modelled by the following 
one-dimensional nonlinear Schrodinger equation [6] 

, BE' d 2 E' l9 n 2 „, , 

Oz' Ox' Uq 

where E' is the envelope of the electric field, k is the wave number in a linear 
medium with n=no, x' and z' are spatial coordinates, and n 2 represents the 
nonlinear part of the refractive index. Within certain power constraints, the 
solution of Eq. ([I]) is a spatial soliton which can propagate in a slab waveguide 
of width much larger than the transverse mode size of the soliton. 

Using the following changes of variables E'=no/kn 2 u, z'=2kt, x'=x, Eq. 
can be written in the following dimensionless form 

iu t = -u xx + Vu, x E V, t > (2) 

where V=— \u\ 2 and T> is the spatial domain. 

Equation g can be written in the standard, dimensional form nonlinear 
Schrodinger equation by transforming the indepedent variables as indicated 
in [3]. In this letter, however, the dimensionless Eq. ^ is used to draw 
an analogy between the propagation of a soliton in a waveguide and the 
motion of a particle in a nonlinear potential wave. This analogy allows to 
determine the dimensionless quantum momentum and quantum energy of 
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and the dimensionless forces on the soliton as follows. The linear quantum 
momentum is defined by the one-dimensional operator [7j 

whose mean value is given by the following time-dependent expression 

(P) = ^ = - M U ^ X , ( 4 ) 

where u* denotes the complex conjugate of it. 
The quantum momentum density is defined as 

p(x,t) = -iu*u x . (5) 

The linear quantum energy is defined by the operator 

whose mean value is 

(E) = = Jf UtdX . (7) 

\u,u) } v u*udx 

while the local energy density is 

e(z,t)=i^-^. (8) 

where g is a factor which has been introduced in order to unify our treat- 
ment with that of Bian and Chan [3] even though these authors used the 
dimensional nonlinear Schrodinger equation to analyze the forces on solitons 
propagating in infinite lines. 
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A physically and mathematically consistent use of the quantum analogy 
can be obtained by employing the following renormalization factor 



Q 



[ u*udx (9) 
Jv 



which coincides with the first invariant of the initial-value problem of the 
nonlinear Schrodinger equation. Since depending on the applications of the 
nonlinear Schrodinger equation and boundary conditions, the first invariant 
may not remain constant, a value g=l is used in this work. 
The energy density can be split as 

e(x, t) = - U Uxx + — = e k (x, t) + e v (x, t) (10) 
Q 

where the kinetic and potential energy densities are, respectively, 

/ u u xx , , \u\ 

e k {x,t) = , e v (x,t) = . 

q g 

From the potential energy density, the following nonlinear force density 
that produces self-focusing on solitons is obtained 

d d /|d 4 \ 

/•(*.') = -£«.<*■*) = S (U-J (ID 

while, from the kinetic one, the following diffraction force density that is 
responsible for the diffraction effect on the soliton is obtained 

The energy and force densities defined above have complex values. In 
quantum mechanics, these complex values do not represent any problem be- 
cause only the mean values of the energy and force can be physically measured 
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and these values are always real due to the Hermitian quantum operators. In 
order to obtain physical insight from the energy and force densities defined 
above, it is convenient to use their real values. 

The force densities may be employed to determine their effects on solitons. 
Bian and Chan [3j used the filter g=u*u in order to assess that a soliton can 
be considered as the result of two opposite phenomena, i.e., the diffraction 
and the self-focusing caused by dispersion and nonlinearity, respectively, for 
the initial- value problem of Eq. (|2]). Note that the force densities obtained 
by Bian and Chan are somewhat artificial and are not consistent with the 
quantum mechanics analogy since their renormalization factor is a function 
of both space and time, and, therefore, affects the values of these densities in 
a different manner depending on the soliton location at each instant of time. 

The numerical results of Bian and Chan [3] for the initial- value problem 
of the nonlinear Schrodinger equation indicate that the total force density, 
i.e., the sum of the diffraction and nonlinear force densities, vanishes for 
the 1-soliton solution. For the general exact iV-soliton solution [HI [9] , their 
definition yields a non-zero total force which is controlled by the nonlinear 
one indicating that solitons are compressed while propagating. 

In this letter, finite planar waveguides in a symmetric, finite interval 
V=[—L,L] subject to the following homogeneous boundary conditions are 
considered 

u(-L,t) + ^u x (-L,t) = 0, u(L,t)+-fu x (L,t) = 0, t>0. (13) 

The values 7=0 and oo correspond to homogeneous Dirichlet and Neumann 
boundary conditions, respectively. The nonlinear Schrodinger equation with 
homogeneous Robin boundary conditions is a skew-symmetric problem in x 
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except for the limiting cases 7=0 and 00 for which it is symmetric; therefore, 
the interaction of a soliton with the left boundary is different from that 
with a right one if 7 is finite and different from zero. In this letter, mixed 
boundary conditions corresponding to 7=1 are considered, and the results for 
the Dirichlet, Neumann and Robin boundary conditions are compared with 
those for the initial-value problem and with those for the following periodic 
boundary conditions 

— (x,t) = —(x+2kL,t), Vn>0, keZ, x e V = [-L, L], t > 0. 

(14) 



3 Presentation of Results 

In this section, some sample results (cf. Figs. 1-5) that illustrate the non- 
linear force density, the real part of the diffraction force density, the total 
force density, and the real part of the momentum density are presented as 
functions of space and time for the four types of boundary conditions con- 
sidered in this letter. Figures 1-5 also show the space-time isocontours of 
the three-dimensional data presented in these figures. The results presented 
in Figs. 1-5 were obtained by means of a Crank-Nicolson finite-difference 
method and correspond to an interval of L=50, amplitude and velocity of 
the soliton equal to one, initial position and initial phase of the soliton equal 
to zero, and spatial and temporal step sizes equal to 0.25 and 0.01, respec- 
tively. These figures show that the nonlinear and diffraction force densities 
have an S'-shape, that the nonlinear force density is larger than the diffrac- 
tion one when the soliton is far away from the boundaries, and that the 
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momentum density has a bell-shape similar to that of the soliton amplitude. 

The boundary does not affect the force and momentum densities on the 
soliton for the periodic boundary conditions as shown in Fig. 1. However, the 
force densities are strongly affected by the boundaries for the Robin boundary 
conditions and the momentum changes sign since the soliton rebounds from 
the boundary, recovering the shape that it had prior to the collision. Figure 2 
corresponds to Dirichlet boundary conditions and shows that the S-shape of 
the force densities is maintained during the collision process although very 
large forces are reached in the collision process since the soliton steepens 
near the boundary because its amplitude is zero there; the momentum at the 
boundary is exactly zero and changes sign smoothly. 

Figure 3 corresponds to homogeneous Neumann boundary conditions and 
illustrates the very large values reached by the force densities and the change 
in their shapes during the collision; the total force density exhibits extrema 
and the diffraction force exceeds the magnitude of the nonlinear one at the 
boundary. 

Since the nonlinear Schrodinger equation with Robin boundary condi- 
tions is not a symmetric problem, the interaction of a soliton with the left 
boundary is expected to be different from that with the right one. This is 
clearly illustrated in Figures 4 and 5 which correspond to the first and second 
collisions, i.e., a collision with the right boundary followed by another one 
with the left boundary. 

Figure 4 indicates that the diffraction force is larger than the nonlinear 
one only near the boundary resulting in a positive total force. The momen- 
tum density behaves similarly to the Neumann case, except that its value 
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at the right boundary is small but different from zero. Figure 5 shows that 
the nonlinear (diffraction) force density exhibits a plateau near the bound- 
ary during the collision before changing its S-shape and reaching a negative 
(positive) extremum at the boundary. Figure 5 indicates that the nonlinear 
force is always greater than the diffraction one when soliton collides with the 
left boundary indicating that the collision with the left boundary is simi- 
lar to that observed with Dirichlet boundary conditions (cf. Fig. 2). The 
isocontours presented in Figures 4 and 5 indicate that the soliton becomes 
closer to the right boundary than to the left one. 

4 Conclusions 

A quantum mechanics analogy is used to determine the forces on and the 
energies of solitons in finite, nonlinear planar waveguides subject to periodic 
and to homogeneous Dirichlet, Neumann and Robin boundary conditions in 
finite intervals. For all the boundary conditions considered in this letter, it 
has been shown that the nonlinear, diffraction and total force densities have 
S'-shape profiles and recover the values that they had prior to the interaction 
with the boundary. The interaction of the soliton with the boundary is char- 
acterized by force densities which increase as the hardness of the boundary 
conditions is increased. 
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Figure 1 Classical (F cl ) and quantum (F qil ) force densities. (Left: infinite waveguide; right: finite waveguide subject to periodic 
boundary conditions) 




Figure 2 Classical (F d ) and quantum (F qu ) force densities. (Left: finite waveguide subject to homogeneous Dirichlet boundary 
conditions; right: finite waveguide subject to homogeneous Neumann boundary conditions) 
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